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ABSTRACT 


A  binary  mixture  theory  with  microstructure  is  constructed  for  unidi regionally  fiber-reinforced 
elastic  composites.  Model  construction  is  based  on  an  asymptotic  scheme  with  multiple  scales  and  the 
application  of  Reissner’s  new  mixed  variational  principle  (1984).  In  order  to  assess  the  accuracy  of  the 
model,  comparison  of  the  mixture  model  predictions  with  available  experimental  data  on  dispersion  of 
harmonic  waves  is  made  for  boron/epoxy  and  tungusten/aluminum  composites.  Formulas  for  the 
effective  moduli  are  also  presented,  and  the  results  are  compared  with  test  data  and  other  available 
predictions. 


1.  Introduction 


With  the  advent  of  high  strength  and  stiffness  fibers  such  as  boron  and  carbon,  and  the  develop¬ 
ment  of  techniques  for  binding  such  materials  to  plastic  or  metal,  fibrous  composites  have  become 
important  elements  of  modern  structures.  Such  composites,  due  to  their  microstructural  heterogeneity, 
may  exhibit  response  phenomena  for  some  environments  that  are  not  observed  for  homogeneous 
materials.  An  example  of  these  phenomena  for  dynamic  environments  is  wave  dispersion,  and  under¬ 
standing  of  which  is  important  both  from  the  standpoints  of  direct  response  prediction  and  indirect  ana¬ 
lyses  associated  with  such  topics  as  nondestructive  testing.  For  fibrous  composites,  wave  dispersion  has 
been  amply  demonstrated  via  ultrasonic  techniques  by  such  investigators  as  Tauchert  and  Guzelse 
(1972),  and  Sutherland  and  Lingle  (1972). 

Simulation  of  response  phenomena  associated  with  the  material  microstructure,  such  as  wave 
dispersion,  requires  a  higher-order  continuum  description.  Several  such  models  have  been  proposed, 
some  phenomenological,  some  nonphenomenological. 

A  higher-order  continuum  model  which  simulates  wave  dispersion  was  first  proposed  by  Achen- 
back  and  Herrmann  (1968)  for  unidirectionally  fiber-reinforced  composites.  This  theory,  called  the 
"effective  stiffness  theory",  has  been  further  studied  and  applied  to  fibrous  composites  by  Bartholomew 
and  Torvick  (1972),  Hlavacek  (1975),  Achenback  (1976),  and  Aboudi  (1981).  The  aforementioned 
work  concerned  linear  materials.  By  modifying  the  original  methodology,  Aboudi  (1982,  1983) 
extended  the  linear  model  to  account  for  inelastic  responses  of  the  composite  constituents. 

In  addition  to  the  effective  stiffness  modeling  concept,  a  mixture  approach  has  been  followed  by  a 
number  of  investigators.  A  phenomenological  version  of  this  model  type  was  adopted  by  Martin,  Bed¬ 
ford  and  Stern  (1971).  Deterministic,  nonphenomenological  mixture  theories  were  introduced  by 
Hegemier,  Gurtman  and  Nayfeh  (1973),  Hegemier  and  Gurtman  (1974),  and  Murakami,  Maewal  and 
Hegemier  (1979).  Although  capable  of  simulating  nonlinear  component  responses  and  interfacial  slip, 
this  work  was  limited  to  waveguide-type  problems.  This  limitation  was  removed  in  the  mixture  theory 
developed  for  laminated  composites  by  Hegemier,  Murakami  and  Maewal  (1979),  and  Murakami, 
Maewal  and  Hegemier  (1982).  In  their  papers,  it  was  demonstrated  that  the  mixture-type  model  was 
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capable  of  simulating  harmonic  wave  dispersion  in  laminated  composites  more  accurately  than  the 
effective  stiffness  theories.  Further,  the  mixture-type  model  requires  fewer  governing  equations.  The 
accuracy  and  efficiency  of  the  mixture  theory  is  due  to  the  use  of  appropriate  displacement  and  stress 
microstructural  fields,  and  a  judicious  smoothing  technique.  These  are  obtained  by  an  asymptotic  pro¬ 
cure  with  multiple  scales.  This  procedure  yields  a  series  of  microboundary  value  problems  (MBVP’s) 
defined  over  a  unit  cell,  which  in  turn  represents  the  (periodic)  microstructure  of  a  composite.  The 
lowest  order  version  of  the  MBVP  method  is  equivalent  to  the  "0(1)  homogenization  theory"  summar¬ 
ized  by  Bensoussan,  Lions,  and  Papanicolaou  (1978),  and  Sanchez-Palencia  (1980).  The  latter,  while  it 
generates  appropriate  static  moduli,  is  nondispersive.  Simulation  of  wave  dispersion  requires  at  least  a 
theory  which  is  classified  as  an  0(e)  homogenization  theory  in  which  e  denotes  the  representative  ratio 
of  micro-to-macrodimensions  of  a  composite. 

To  date  an  0(e)  mixture  theory  has  not  been  constructed  for  fibrous  composites  subject  to  arbi¬ 
trary  wave  motion.  Construction  and  validation  of  such  a  3D  model  for  unidirectional  binary  compo¬ 
sites  with  periodic  microstructure  are  the  objective  of  this  paper.  To  facilitate  this  task,  the  asymptotic 
procedure  with  multiple  scales  noted  previously  is  combined  with  a  variational  technique  (Murakami, 
1985).  Following  development  of  the  basic  equations,  the  dispersion  of  time-harmonic  waves  is  studied 
and  the  results  are  compared  with  experimental  data  for  boron/epoxy  (Tauchert  and  Guzelse,  1972) 
and  tungusten/aluminum  (Sutherland  and  Lingle,  1972)  composites.  The  good  correlation  obtained 
with  experimental  data  indicates  that  the  proposed  mixture  model  furnishes  a  basic  tool  by  which 
dynamic  responses  of  elastic  composites  can  be  investigated.  While  the  model  construction  procedure 
is  applicable  to  inelastic  component  response  and  interface  slip,  extension  and  investigation  of  the  non¬ 
linear  problem  is  deferred  to  later  publications. 


2.  Formulation 

Consider  a  domain  V  which  contains  a  uniaxial  periodic  array  of  fibers  embedded  in  the  matrix,  as 
shown  in  Fig.  1.  Let  a  rectangular  reference  system  3c j,  x2, 3c j  be  selected  with  x{  in  the  axial  direction 
of  the  fibers.  In  the  x2 ,  jfj -plane,  a  typical  cell  that  represents  the  geometrical  microstructure  of  the 
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composite  is  shown  in  Fig.  2  for  a  hexogonal  array. 

For  notationa!  convenience  forms  (  )fa),  a  -  1,2  denote  quantities  associated  with  material  a 
with  a  «  1  representing  fiber  and  a  -  2  matrix.  Cartesian  indicia!  notation  will  be  employed  in  which 
Latin  indices  range  from  1  to  3  and  repeated  indices  imply  the  summation  convention  unless  otherwise 
stated.  In  addition,  the  notations  C)^  =  8(-)/0x,  and  C)j  =  dO/dT  will  be  employed  in  which  7 
represents  time.  Quantities  of  the  form  (“)  and  (  )  denote  dimensional  and  nondimensional  variables, 
respectively. 

The  governing  relations  for  the  displacement  vector  ul<a)  and  the  stress  tensor  in  the  two 
constituents  are: 

(a)  Equations  of  motion 

— (<*)  /1\ 
“ijt  >  ™  Oy  VU 

where  p(a>  is  the  mass  density; 

(b)  Constitutive  relations 

-  Afa)8„-  +  (2) 

where  X(<*\  /I*"1  are  Lame’s  constants,  elja)  is  the  infinitesimal  Cauchy  strain,  and  80  is  the  Kronecker 

delta; 


(c)  Strain-displacement  relations 


^  -  4  Gtf  ’  +  > 


(d)  Interface  continuity  relations 


V*  —  u,<2>  ,  —  of’*;0  on  (4) 

where  i»{!)  =  0  on  the  fiber-matrix  interface  &  ; 

(e)  Initial  conditions  at  7  —  0  and  appropriate  boundary  data  along  the  boundary  d  V. 

Conditions  (a)  -  (e)  define  a  well  posed  initial  boundary  value  problem.  However,  due  to  the 
large  number  of  fiber-matrix  interfaces  the  direct  solution  to  this  problem  is  extremely  difficult.  The 
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objective  of  the  subsequent  analysis  is  to  alleviate  such  difficulties  by  deriving  a  set  of  partial  differential 
equations  with  constant  coefficients  whose  solution  can  be  utilized  to  approximate  the  solution  of  the 
problem.  To  this  end,  it  will  be  convenient  to  nondimensionalize  the  basic  equations  by  using  the  fol¬ 
lowing  quantities: 

A  typical  macrosignal  wavelength 

A  typical  fiber  spacing  or  cell  dimension 


C(m)  i  P(m ) 

^{m)  =  P(m)  C(m) 
^(m)  —  A/C(m) 

«  =  A/A 


reference  wave  velocity  and  macrodensity 
reference  modulus 
typical  macrosignal  travel  time 
ratio  of  micro-to-macrodimensions. 


With  the  aid  of  the  above  notation,  nondimensional  variables  are  now  introduced  according  to 


(*1  ,  x2  ,  x3)  -  Cx,  ,  x2  ,  x3)/\  ,  t  -  , 


(X^- (X£)(*y£(m)  ,  pfa)*p(a)/p(w)  .  (5) 


With  the  variables  defined  according  to  (5),  the  material  properties  are  seen  to  be  periodic  in  the 
x2,  xj-plane  in  which  the  periodicity  of  the  fiber  lattice  structure  may  be  defined  by  the  cell.  It  is 
expected  that  stress  and  deformation  fields  will  vary  significantly  with  respect  to  two  basic  length  scales: 
(I)  a  "globaP  or  "macro"  length  typical  of  the  body  size  or  loading  condition,  and  (2)  a  "micro"  length 
typical  of  "cell"  planar  dimensions.  Further,  it  is  expected  that  these  scales  will  differ  by  at  least  one 
order  of  magnitude  in  most  cases.  This  suggests  the  use  of  multivariable  asymptotic  techniques  (Ben- 
soussan,  Lion  and  Papanicolaou,  1978,  Hegemier,  Murakami  and  Maewal,  1979,  Sanchez-Palencia, 
1980).  This  approach  commences  by  introducing  new  independent  microvariables  according  to 

■xi*  =  x,lt  .  (6) 

Therefore,  all  field  variables  are  considered  to  be  functions  of  the  microvariables  x\  and  xj,  as  well  as 
the  macrovariables  x, ,  i  -  1-3: 

/(*!,  X2,xh  f)-/'(x|,  X2,  *3,  x\,x\,  t  ;  e)  (7a) 
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Spacial  derivatives  of  a  function  /  then  takes  the  form 


f  (x*  ,  r)  -  f  (x* ,  xj,  t  ; «) 


+  «  dx;r  °Ck'x;,t  'e) 


(7b) 


where  d(  )/9xJ  =  0.  By  introducing  the  notation  (  ),.  =  d(  )/dx *  equation  (7b)  can  be  rewritten 
as: 


/j-^  +  7^- 


(7c) 


In  the  sequel  /*  will  be  written  as  /  for  notational  simplicity. 

The  operations  (7),  when  applied  to  all  field  variables,  lead  to  the  following  "synthesized"  govern¬ 
ing  field  relations: 

(a)  Equations  of  motion 


o-fp)  +  1  _.(«).•(»> 

"jij  T  „  vjij  P 


•  to) 


.(a) 


(8) 


(b)  Constitutive  relations 

< r P  -  Xfa)  ; 


(c)  Strain-displacement  relations 


».(*) 


^’+7  tyy+ufti)  ; 


(9) 


(10) 


(d)  Interface  continuity  conditions 

,  ajV’  “  on  <f  .  (11) 

At  this  point,  the  variation  of  field  variables  which  satisfy  the  periodicity  with  respect  to  x*  is  assumed. 
According  to  this  condition  field  variables  take  equal  values  on  opposite  sides  of  the  cell  boundary.  The 
premise  allows  one  to  analyze  a  single  cell  in  an  effort  to  determine  the  distribution  of  any  field  variable 
with  respect  to  the  microcoordinates  x*.  The  x*  -periodicity  condition  is  motivated  by  the  Floquet  and 


Block  theorems  (Brillouin,  1946)  for  harmonic  wave  in  periodic  structures.  Certainly,  it  eliminates 
boundary  layer  effects.  However,  it  is  expected  to  provide  a  good  model  for  the  global  wave 
phenomena  in  fibrous  composites  with  periodic  microstructure. 


For  the  construction  of  a  mixture  model  it  is  convenient  to  cast  the  field  equations  in  a  variational 
form  by  using  the  Reissner  new  mixed  variational  principle  (Reissner,  1984).  In  the  Reissner  varia¬ 
tional  principle  the  variations  of  displacement,  strain  with  (10)  as  definition  and  transverse  stresses,  i.e., 
all  stress-components  except  cr/f*,  are  considered.  Thus,  it  is  convenient  to  rewrite  the  constitutive 
relation  (9)  in  terms  of  the  axial  strain  e}f)  and  the  transverse  stresses: 


orft*  =  (X  4  2fjL)(a)e[ f>  4  Xfa) 


(•••)  +  effl  (•••)}  , 


e#  <  • 

•  ) 

( • 

•  ) 

(X  4  3*)(a)  Xfa) 

Xfa)  (X  +  $i)(o) 


-l 


0"22| 

CT33] 


(a) 


X<a)e,<f> 


[2eft>  ( •  •  •  )  ,  2etf  ( •  •  •  )  ,  left >  ( •  •  •  )] 


_  1 


4)  Itrjy’  ,  crjf’  , 


(12) 


Using  the  equations  of  motion  (8),  Gauss’  theorem,  and  the  x' -periodicity  condition,  it  can  be 
demonstrated  that  the  Reissner  mixed  variational  principle,  applied  to  the  synthesized  fields  by  the  mul¬ 
tivariable  representation,  takes  the  form: 


4  28e3(?)  o-jf*  4  28e,(?) 


4  4  -  ufol  -  e2<5>  (•••))  4  4-  -  u$.  -  eft  (•••)) 


4  85-  j?’(w $  4  u$  +  j  wft'  +  j  ~  2e2(J)  (•••)) 
4  4  uft  4  j  ufo>.  -  2ejf)  (•••)) 


k.'-  .  ,n  : 


vVW v^v*. V.V*V^ V«\ * •.**-  v  > V.*V* VA v\v V*V' *1' 
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+  Sd'ft’fe  fe*  +  "ft  +  j  u  fjl  -2e{^  (■■•)) 
-  +J  |(8«4(2)  -  &u,(l))Ti  +  87^(u,(2)  —  t4(1))J</s* 

UH  HJ 


dx\  dx\ 


dx\  dxi  cb C3 


8u,<a)(—  p(a)Ui[°})dx j  cb cH  dxj  dx2dx} 


+  ff  <  I  ff  Sw/0’  7?a)  dx\  dx]  )dA  , 
a‘‘  A& 


(13) 


where  A (o)  denotes  the  xj ,  x 3 -domain  of  the  cell  occupied  by  material  a  (Fig.  2),  <x,‘Q)  is  used  for  the 

V 

approximate  transverse  stresses,  Tfa)  denotes  the  traction  vector  on  the  surface  9  VT  where  the  traction 
is  specified,  ds'  is  an  infinitesimal  line  element  on  & ,  and  dA  is  an  infinitesimal  surface  element  on 
the  boundary  of  V  :3  V.  In  (13)  basic  variables  are  the  displacement  u,<a),  the  transverse  stresses 

V 

and  the  interface  traction  vector  T".  The  Euler-Lagrange  equations  of  (13)  include  (8a),  (11a),  (12), 
and 


on  &  . 


(14) 


The  above  variational  equation  (13)  furnishes  a  tool  with  which  a  mixture  model  can  be  obtained 
with  appropriate  trial  displacement  and  transverse  stress  fields.  The  basic  requirement  for  the  variables 
is  the  x* -periodicity  condition  on  the  cell  boundary  dA  .  The  microstructural  variation  of  the  trial  func¬ 
tions  can  be  obtained  by  the  asymptotic  procedure  (Murakami,  Maewal  and  Hegemier,  1981). 


3.  Asymptotic  Analysis 

The  premise  that  the  composite  macrodimension  is  much  larger  than  the  microdimension, 
«  «  1,  and  the  form  of  scaled  equations  (8)  and  (10),  suggest  the  expansion  of  the  dependent  vari¬ 
ables  in  the  asymptotic  series: 


U/TyJ^O^X,*,  t;«)=  I  €''{w((,,,«Ty(n,)(a)(x*,X/’,  r) 

n—o 


(15) 
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If  (15)  is  substituted  into  (8)-(ll)  and  the  coefficients  of  different  powers  of  e  are  equated  to  zero,  a 
sequence  of  problems  defined  on  the  cell  is  obtained.  The  first  of  the  equations  in  this  sequence  fur¬ 


nishes 


,  vjftUy  -  0  .  (16) 

Equation  (16a)  implies  that  u,\a0]  is  independent  of  xj  and  yields  with  the  zero-th  order  expansion  of 


(11a): 


”  Vl{9)  (x*  ,  t)  . 


The  remaining  systems  of  equations  obtained  from  (8)- ( 10)  are,  for  n  ^  0: 


<Tji\>n),j  '  ** ji\n  )  “  ij\n  )  » 


-  4  + + + uftiw)  ■ 


To  be  added  to  the  foregoing  are  the  interface  conditions  and  the  x'  -periodicity  conditions  for  n  ^  0: 

“  “iff)  ,  <rX)p}l)  =  <r}ft„)v}l)  on  &  ,  (21) 

utfj)  and  <rj$,)Pj2)  are  x*-periodic  on  dA  .  (22) 

The  first  set  of  microboundary  value  problems  (MBVP's)  for  and  u,^},  called  the  0(1) 
MBVP’s,  is  defined  by  (16b),  (18b),  (19)-(20),  (21b),  (22b)  with  n  =  0,  and  (21a),  (22a)  with  n  =  1. 
The  0(1)  MBVP’s  are  excited  by  U^0)j.  Similarly,  a  sequency  of  MBVP’s  is  defined  for  each  n  from 
(18M22).  With  appropriate  integrability  and  normalization  conditions,  higher  order  terms  may  be 
computed  by  solving  the  MBVP's.  In  particular,  the  0(1)  MBVP’s  are  the  ones  solved  for  the  0(1) 
homogenization  theory  proposed  by  Bensoussan,  Lion  and  Papanicolaou  (1978)  and  Sanchez-Palencia 
(1980),  and,  also,  form  the  basis  of  the  mixture  theory  which  may  be  classified  as  an  0(e)  homogeniza¬ 
tion  theory.  The  asymptotic  approach  yields  the  microstructures  of  displacement  and  stress  fields  after 
solving  a  multitude  of  MBVP's  which  are  complicated. 
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In  order  to  use  the  approximate  solutions  of  the  MBVP's  in  the  course  of  developing  a  mixture 
model,  and  to  ease  the  burden  of  solving  the  MBVP's  exactly,  a  variational  procedure  was  adopted  by 
Murakami  (1985)  for  laminated  composites  with  the  help  of  the  Reissner  new  mixed  variational  princi¬ 
ple  (Reissner,  1984).  A  similar  approach  is  adopted  here  for  fibrous  composites.  To  obtain  the  lowest 
order  mixture  theory  by  using  (13),  it  is  necessary  to  obtain  trial  displacement  and  transverse  stresses 
to  0(e).  In  the  sequel,  the  trial  functions  are  obtained  for  a  hexagonal  cell  with  a  concentric  cylinders 
approximation  as  shown  in  Fig.  2.  In  Fig.  2,  (r  ,9)  are  micropolar  coordinates: 

r  —  y/x\2  +  xj2  ,  tan0«x^/xj  »  (23) 

by  which  r  -  1  constitutes  the  cell  boundary  and  r  —  y/n^\  denotes  the  interface  & .  The  quantities 
/i(“)  indicate  the  volume  fraction  of  material  a  and  satisfy 

+  „«>  -  1  .  (24) 

In  terms  of  the  polar  coordinates  the  x  ‘-periodicity  conditions  for  a  hexagonal  cell  with  the  concentric 
cylinders  approximation  reduce  to  the  form: 

f(xrft,t)  - /Gi  ,r,ir  +  9,t)  at  r  -  1  .  (25) 


4.  Trial  Displacements  and  Transverse  Stresses 

The  0(1)  stress  and  0(c)  displacement  fields  are  obtained  by  solving  the  0(1)  MBVP’s  which  are 
defined  by  (16b),  (18b),  (19)-(22)  and  (24).  These  MBVP's  are  excited  by  Uno)j.  The  exact  solution 
of  is  furnished  in  the  Appendix.  For  the  mixture  formulation  it  is  convenient  to  introduce  an 
0(c)  displacement  variable  which  represents  Un0)j  +  (Jn0)j  according  to: 

Sibck,t)  =  J  u,(a)vjX)ds'  -  -7/  (26) 

where  ^(-rr)  is  the  area  of  the  cell.  Due  to  the  fact  that  uftf  is  excited  by  U^o)j  +  £/,<<,),  one 
obtains 

3  2 

. 


(27) 


-  12  • 


i m 


Equation  (27)  can  also  be  obtained  if  one  substitutes  the  exact  in  the  Appendix  into  (26)  and 
eliminates  U,{o)J.  To  render  the  analysis  tractable,  it  is  preferable  to  utilize  an  approximate  form  of 
the  exact  solution  for  The  exact  solution  indicates  that  the  following  form  of  the  0(e)  displace¬ 
ment  yields  a  good  approximation: 

*4^1  r Xj.  ,r)  -  ^  (xk  ( r )  cos  9  +  S,  Cx*  ,f)gfa)  (r)  sin  9  (28a) 

where 

s(,)  '  g(2)  (r) "  -Jir  (~ r  +  7}  •  (28b) 

Anticipating  the  0(e2)  difference  of  the  average  of  up  on  A (a),  equations  (17)  and  (27)  yield  the  fol¬ 
lowing  trial  displacement  field: 

uP  (x*  ,Xj. ,/)  -  UP  Cx*  ,r)  +  e  (x*  ,/)  (29) 

where  is  defined  by  (28).  Equations  (29)  and  (28)  indicate  that  the  mixture  displacement  vari¬ 
ables  are  Up,  U,{2),  S,  and  S,  with  the  constraint  (27). 

By  using  (29)  in  (19)  with  n  —  0  and  considering  the  0(e2)  differences  of  the  average  transverse 
stresses,  the  0(1)  trial  stress  field  may  be  expressed  as: 
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In  order  to  define  the  0(e)  trial  stress  field  it  is  convenient  to  define  the  0(e)  stress  variable 
according  to 

P,  w  ,«>  =  ±  ±  S  ■  (31) 

If  one  integrates  (8a)  over  /4<o)  and  utilizes  the  x'  -periodicity  condition,  one  obtains  the  mixture 
momentum  equations: 


mm 


(32) 


n 40  +  (-1  )“+1  /»  -  n0*’  p(“)  i^o) 

where  the  average  operation  is  defined  by 


f(aa)0ckj)  =  dx]dx\  (33) 

From  (32)  it  can  be  seen  that  P,  represents  an  interaction  body  force  between  the  two  constituents 
across  the  interface.  Also,  the  form  of  (32)  with  P,  defined  by  (31)  satisfies  the  integrability  condition 
adopted  by  the  0(1)  MBVP's  for  which  are  defined  by  (18)-(20)  with  appropriate  n’s  and  (31). 
As  an  0(e)  trial  stress  field  which  satisfies  (31b)  one  may  use  the  following  approximate  fields: 


^22(1) 

<o) 

.  1 
”  4 

3 

ll 

O' 33(1) 

,t)gla)  ( r )  cose 

1 

+  Pi(xkj)gia)(r)  sin© 

3 

Or23(l) 

1 

1 

fa) 

o^Kl) 

-  \  Pi(xkj)g(aHr) 

sin  e 

O' 12(1) 

cose 

(34a) 


(34b) 


As  a  result,  the  trial  transverse  stresses  are  expressed  as: 

=  au{o)(xk  ,X/V)  +  e  6v/(i)k*  r*/V)  (35) 

where  afy ,  and  o',]"})  are  defined  by  (30)  and  (34),  respectively. 


5.  Mixture  Equations 

By  substituting  the  displacement  and  transverse  stress  trial  functions  defined  by  (29)  and  (35), 
respectively,  into  the  Reissner  variational  equation  (13),  one  obtains  the  following  relations  as  the 
Euler-Lagrange  equations: 

(a)  Equations  of  motion 

+  (-1  )“+1^  -  ,  /  *  1-3  ,  (36) 

Mj,j+ -^(<ri2o)-<rUa)+ Rtf')- l  1,2  , 


(37a, b) 


(37c,d) 


Mju  +  ji  ^  ~  W/0>  +  *&)  -  /  Si, a  ,  /  -  1,3  , 
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«-l  4  4«(2)  n(2) 


(b)  Constitutive  relations 
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M22-hP2,M33- 4  *2.^12- 4  4  **J  - 


2  *  *  ’  w  2  1  ’  ^  2 

2  2  2  3  3 


M33  —  hP3  ,  M3\  —  -j  Pi  ,  M23  —  A/23  “  A/31  —  A/23  —  A/12  —  0 


(43) 


where  it  is  understood  that 


2  2  3  3 

My  ™  A/,,  ,  A/^  *  Mjj 


(44) 


*2(P  = 


r8,//i‘2)  ;  Rtf  =  (ftf/2  +  ttf  >/n(2)  ,  Pi?'  =  -  ttf/n(2)  , 
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(45) 


and 


V«(2)  , 


-  (X  +  M)(»(Sj  ~  S3))/«<2)  ,  'if  “  ~  M(2)«2  +  53)/"(2) 


ttf  -  -  (x  +  /u.)<2)(2  s2)/» 
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(46) 


The  remaining  constitutive  relations  associated  with  o-|(f)  are  obtained  from  (12a);  the  results  are 


<rfta)  -  (X  +  2/a)<“)tfft)  +  Xfa) 


£/&>  +  C/fc’  +  (- 1 )“+ 1  (12  +  S3)/n  <“  >  , 


(47) 
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The  associated  boundary  conditions  are  on  d  V 

n'a)<rfca)vj-T}ap)  .  or  SC/J<“>  —  0 


1- 1,2,3 


(49) 


2  2k 
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(Mj2  +  Mjj)vj  ■  Tj+  7j  or  8  Sj  ■  0  .  (50c) 

where 

T!*''=±f£Tfa)dx’2dxl  . 

€  r(  ,  rj  -  -7  I  ,  Sine)  dx\  dx\  .  (51) 

A  o-l 

Equations  (36)-(50)  and  the  initial  condition 

Up*  ,  UP  ,  i  ,  S(  ,  at  r  -  0  on  K  (52) 

define  a  well  posed  initial  boundary  value  problem  with  respect  time  t  and  the  macrocoordinates  xk . 


6.  Harmonic  Wave  Dispersion  Spectra 

In  an  attempt  to  test  the  accuracy  of  the  mixture  model,  the  phase  velocity  and  group  velocity 
spectra  of  the  mixture  theory  have  been  compared  with  available  experimental  data  for  time  harmonic 
waves.  For  the  comparison  harmonic  waves  which  are  propagating  at  an  arbitrary  angle  of  incidence  in 
a  full  space  of  the  following  form  are  considered: 

,  U{2)  ,  UP  ,  U{2)  ,  UP  ,  UP  ,  Sy/ik  ,  SJtk  ,  ISJik  ,  SjikV 

—  exp  | ik  (x\  cos  #  +  x2  sin  #  cos  9  +  x3  sin  #  sin  9)  —  i'o><J  If  (53) 

where 


If  -  [i/P  ,  UP  ,  Up  ,  i/P  ,  UP  -  UP  J|  i  J,  ,  ljlr 


(54) 


and  {  lr  denotes  the  transpose  of  [  ].  In  (53)  UP  and  s,  are  constant  amplitudes,  k  denotes  the 
wave  number,  w  represents  angular  frequency,  #  is  the  azimuth  measured  from  the  X\  axis,  and  9  is 
the  longitude;  the  direction  of  the  wave  propagation  may  be  best  represented  by  the  wave  vector  k: 


/r—  k (cos #  ,  sin#  cos 9  ,  sin#  sinfl)r 


(55) 


Substitution  of  (S3)  into  (36)  and  (37),  which  are  expressed  by  the  displacement  variables  with 
(40). (46),  yields  an  eigenvalue  problem  for«<u  of  the  form: 

[*]£/-  («o >)*[MW  (56) 

where  [A"  ]  and  [M]  and  are  11  x  11  real  symmetric  matrices,  the  elements  of  which  are  functions  of 
the  mixture  constants  and  the  wave  vector.  Furthermore,  [M]  is  a  diagonal  matrix.  Upon  calculation 
of  the  eigenvalue  eat  for  a  given  ek,  one  obtains  the  phase  velocity  Cp  as 

Cp  -  (eo>)/(ek)  .  (57) 


For  each  computed  eigen  pairs  (eto  ,  0k,  k  -  1,2,  •  •  •  11  the  group  velocity 


can  be  obtained  by  taking  the  derivative  of  (56)  with  respect  to  ek: 

+  (cg>)2[A/']  Uk  . 

k 

For  the  kth  eigenpair  equation  (59)  yields 


ir]*#- 


2(e<o)CglM] 


(58) 


(59) 


(«o;)2[A/']}  U 

(Cf  )*  —  - - : - : -  •  (60) 

2(e<v)*(f/7[A/)C0* 

In  the  subsequent  simulation  a  typical  cell  dimension  A  was  chosen  to  be  a  cell  radius  by  introducing 
the  concentric  cylinders  approximation  of  the  equal  area.  The  reference  elastic  modulus  and  density 
used  for  the  scaling  are 

£(m)-  I  «<“>£<“>  ,  p(m)-  r  (61) 

a— 1  o-l 

where  £ <o)  is  Young’s  modulus.  The  dimensional  frequency  v  (Hz )  can  be  computed  from  «o*  by 
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(««  K/£(m/P(m^  (2 IT  A  ) 


(62) 
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Numerical  results  are  presented  for  a  boron-epoxy  composite,  for  which  experimental  results  were 
presented  by  Tauchert  and  Guzelse  (1972)  for  a  waveguide  case  4>  —  0°  and  a  waveflect  case  <f>  «  90°. 
The  material  properties  are  summarized  in  Table  1  in  which  the  values  for  Poisson’s  ratio  are 
estimated.  In  the  simulation  A  was  computed  from  the  fiber  diameter  (-  2n/«(1,A)  which  was 
1.016  x  10*4  m.  The  group  velocity  spectra  for  a  waveguide  case  <t>  =  9  —  0°  are  shown  in  Fig.  3  for 
two  acoustic  modes:  a  "gross"  longitudinal  mode  and  a  "gross"  shear  mode.  In  the  figure  the  same  sym¬ 
bols  as  the  reference  of  Tauchert  and  Guzelse  are  used  for  the  experimental  data  points.  It  is  noted 
that  reasonable  agreement  is  achieved  for  the  waveguide  case  in  which  pronounced  dispersion  is 
observed.  The  group  velocity  spectra  for  a  wavereflect  case  <6  —  90°,  9  —  0°  in  which  the  wave  vector  is 
normal  to  the  fiber  axis  are  shown  in  Fig.  4  with  the  experimental  data.  The  figure  includes  three 
acoustic  modes:  a  "gross"  longitudinal  wave  (P-mode),  a  "grostf’  vertically  polarized  shear  wave  (SV- 
mode),  and  a  "gros^  horizontally  polarized  shear  wave  (SH-mode).  The  sets  of  experimental  data 
correspond  to  the  "gross?  P-mode  and  the  "gross?  SH-mode.  It  is  noted  that  there  are  significant  devia¬ 
tions  from  the  "grostf*  SH-mode,  but  the  overall  agreement  is  not  unsatisfactory  if  one  admits  the  scar¬ 
city  of  the  experimental  data  and  the  difficulties  associated  with  the  measurement  of  shear  wave  veloci¬ 
ties.  It  was  reported  by  Tauchert  and  Guzelse  (1972)  that  a  shear  wave  exhibited  extremely  high 
damping  of  the  pulse.  A  similar  observation  and  the  scatter  of  shear  wave  data  were  reported  by  Sachse 
(1974)  who  conducted  modulus  measurements  of  boron/epoxy  composites  by  using  pulse-echo  tech¬ 
niques.  He  concluded  that  "the  measurements  of  the  present  investigation  indicate  that  shear  waves 
propagating  along  and  across  fibers  in  the  composite  materials  tested  do  not  always  propagate  at  the 
same  speed." 

Sutherland  and  Lingle  (1972)  reported  phase  velocity  measurements  for  tungusten/aluminum 
composites  whose  material  properties  are  shown  in  Table  2.  The  equivalent  cell  radius  A  was  computed 
from  the  given  fiber  spacings  which  yield  the  area  of  a  typical  cell  A  (—  tt  A2)  0.579  x  1(T6  m2.  Fig¬ 
ure  5  shows  the  phase  velocity  vs.  frequency  relation  for  the  "gross"  longitudinal  mode.  A  reasonable 
agreement  is  observed  between  the  experimental  data  and  the  theoretical  prediction. 


7.  Effective  Moduli 

The  0(1)  homogenization  theory  which  yields  the  effective  moduli  of  composites  can  be  obtained 
by  taking  the  limit  of  «  — *  0  and  equating  the  constituents'  displacements 

t/,(1)  -  Uj  (63) 

By  introducing  the  above  constraints,  equations  (36)  yield 


where 


Equations  (37)  yield 


rr  (n)  —  n  (m )  // 

a jij  P  Vi,, 


<Tym)  «  I  n(a)cr<afl)  ,  p(m)-  I 


,  /- 1,2,3 
,  /*1,3  . 


By  eliminating  S,  by  (66),  equations  (65a),  (40)  and  (47)  with  (63)  furnish 


r(m)  _  J^OntJ^On) 


where 


<r(m)  *  fcrfr*  ,  <r2(2m)  ,  o-jf  >  ,  ,  a3(r)  ,  o-,(?)]7'  , 

?<m)  “  ((/|,i  i  Ui,2  *  Uu  ,  f/2,3+  f/3,2  ,  £/31  +  t/i  3  ,  t/1<2  +  i/2IF  ,  (68) 

and  is  the  effective  modulus  matrix  with  transverse  isotropy  due  to  the  concentric  cylinders 

approximation  and  is  defined  in  the  Appendix. 

The  formulas  for  the  effective  moduli  (B2)  are  assessed  by  comparing  the  results  with  the  experi¬ 
mental  data  reported  by  Datta  and  Ledbetter  (1983)  for  boron/aluminum  composites.  The  results  are 
shown  in  Table  3  in  which  the  moduli  computed  from  the  effective  stiffness  theories  for  the  square  cell 
by  Achenback  (1976)  and  for  the  hexagonal  cell  by  Hlavacek  (1976)  are  included  by  using  the  formulas 
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reported  by  Datta  and  Ledbetter  (1983).  The  comparison  has  revealed  that  all  high-order  theories  yield 
almost  similar  results.  It  can  be  easily  shown  that  the  formulas  for  the  effective  moduli  yield  values 
which  fall  between  the  upper  and  the  lower  bounds  obtained  by  Hashin  and  Rosen  (1964)  for  fiber- 
reinforced  composites. 


8.  Concluding  Remarks 

An  asymptotic  mixture  theory  with  multiple  scales  was  applied  to  unidirectionally  fiber-reinforced 
elastic  composites  with  periodic  microstructure.  In  the  model  construction,  Reissner’s  new  mixed  vari¬ 
ational  principle  was  applied  to  the  synthesized  fields  with  multivariable  field  representations.  In  order 
to  assess  the  accuracy  of  the  model  the  mixture  dispersion  spectra  were  compared  with  the  experimen¬ 
tal  data  obtained  for  the  boron/epoxy  composite  by  Tauchert  and  Guzelse  (1972)  and  for  the 
tungusten/aluminum  composite  by  Sutherland  and  Lingle  (1972). 

A  satisfactory  correlation  with  the  experimental  data  indicates  that  the  proposed  mixture  model 
furnishes  a  basic  tool  by  which  dynamic  responses  of  the  composite  structures  can  be  investigated. 
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Appendix 


A.  Exact  of  the  0(1)  MBVP’s: 

“1*71')  *■  C  ^<o)(r){  (£/ KO) 2  +  U 2(0).  1 )  cos  6  +  (£/3(o>,!  +  ^ uo),3^s*n  #}  » 

u$U  “  b0(U 2(0, ,2  +  t/3(0).3  +  A  L/ K0).l  )s(a>  (^)cos  0 

+  ^“’[{gk’OOcosfl  +  k ^ 1  gi>a)(r)  cos  36}(f/2<0>,2-  UM0)j) 

+  (£W  (r) sin  6  +  *<“'  &}a)  (r)sin  36}  W2m, 3  +  i^tou)] 

+  fc2(a)|{3(l  -  K<a))^/a>(r)cos  9  +  (1  +  kla))giff)(r) cos  36}  ((/2(0).2 - 
+  {3(1 -K(“))ft(a)(r)sin6+  (1  +  (r)sin  36} (£/2(ou  +  t/3(0 

03%  “  M^2<0),2  +  ^3(0)3  +  *  ^K0).i)f<“)('-)sin  6 

+  a^)[u<“>(r)sin  6  -  *(a),&}a>sin  36} (-  i/2(0>,2  +  Umj) 

+  {f^OOcos  6  -  k^gfrKr) cos  36}  (tf2(0W  +.  £/3(o>.2>] 

+  ufa)  (/,)cos  6  -  K(“)ft}a>  (r)cos  36}  (t/2(0), 3  +  £/3(0),2)] 

+  62(o)|{3(1  -  K(a))g/a)(r)sin  6  -  (1  +  Kfa))f/}?)(r)sin  36)  (-  £/2(o,,2 

+  {3(1  -  K<a))g,(aHr)cos9-  (1  +  K<ar) )giff}r )  cos  36}  (£/2(0).3  +  ^3<o: 

where 

60  =  {(X  +  /k)(2)  —  (X  4-  pi)(1)}/(2c/i )  ,</]=£  (X+^)<“)/«<,,)  +  nl2) / (n{i) n(2) ) 

ar-  1 

(7  =  —  (/t ( n  —  /x {2))/i/ 3  ,  d3=  L  n<a)/n(a)  +  n(2)/ (ni])n(7))  , 

a-  I 

X  =  (xu>  -  \(2>)/{(A+M)(1>  -  (X+m)(2,J 

,  Kfa)=  (1— K(o>)/(l+K  <o)) 


(X+^l)(o)/Mto) 


giiu(r)  -  gfflir)  -  r3  ,  g^Hr)  -  0 

gii2Hr)  -  -  r}+  r~l  ,  g/fHr)  -(-/•"'  +  r~3)/nl2)  , 

gfiPb)  —  r3  —  (1+k(2>)-2  {4(1— k  <2)+*c  <2>2)r-3  —  3(1— #c<2))2r-1}  ,  (A3) 

In  Eqs.  (Al)  a2a)  and  Ma>  are  obtained  by  solving  the  linear  equations  for  (a j 1 1 ,  bP ,  aj2), 

bp  ]r: 

U)  *  -  B  (A4) 

4x4 

4x1  4x1 

where 

An~  l,A2l-0,A3l-pw/n<n  ,  A4l  -  0 

Au- 3(l-K(»)n(1)J  ,  ^ -  (1+k<»)*<»2 

An“  ~  AA2—  3n(1)/x(1)(l— k(i))  ,  >4  13—  —  1  ,  A  23—  —  k<a)/n<a)  , 

,433-MX2,/«(2){1-k<“,/«(1))  ,  ^43-3K(2V«)/n(,)2 

/414--3n(2)(l  +  n(1)Ml-K(2)) 

An  —  (1  +  k<2))/7(1)2  —  (4  —  3/t(1))(l  —  k(2>  +  k(2)2)/{/j(1)(1  +  k(2))} 

^34—  3/4<2)(l  —  K<2)Mn(l)—  k(2)/«(1)) 

^44*  ~  3/i<2)(l  —  k(2))(/?(1)  +  3Sr(2)/n(1))  +  12^t(2)(l  -  k<2)  +  k<2)2)/{/j0)2(1  +  k<2)>)  ,  (A5) 
and 

Bx  «  B2  -  fi4  -  0  ,  *3  -  -  <m(1)  -  fi(2))/2  .  (A6) 


It  is  interesting  to  note  that  for  most  of  practical  composites  b-f} ,  a  -  1,2  are  small  compared  with 


B.  The  Definition  of  tE(ml]  in  (67) 


£ ii  fu  £"12  0  0  0 

£12  ^22  ^23  0  0  0 

£  12  £33  E  22  0  0  0 
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0  0  0  0  Ess  0 
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where 


Eft'  -  1  n<*'(\  +  V)(a)  -  (x(1)  -  \{1)Yldx  , 


£,(2m)-  I  n(a)X(o)-  <X(1)-X(2)){(X  +  M)(1)-  <\+n)(2))/d 
0—1 
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a— 1 


£i7’  -  (Eft'  -  Eft')/ 2  , 


Eft'  -  I  /!<"V<o)-  k(,)-M(2W 3  , 
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and  where 


d2  -  I  M(“)/«(o)  +  (X  +  M)(2)/(2n<,)/j(2)) 
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Table  1.  Material  Properties  of  the  Boron/Epoxy  Composite  Tested  by  Tauchert  and  Guzelse  (1972) 


Volume 

Fraction  n<a) 

Young's  Modulus 

£<a) 

Poisson's  Ratio 

Mass  Density 

Pfa) 

^Boron 

0.54 

379.2  GPa 

(55  x  ltf  psi) 

0.18 

2682  kg/m3 

(251  x  ir6  lb  sec*/in*) 

(2)Epoxy 

0.46 

5.033  GPa 

(0.73  x  1(P  psi) 

0.40 

1261  kg/m3 

(118  x  10*6  lb  sec 2/in*) 

Table  2.  Material  Properties  of  the  Tungsten/ Aluminum  Composite  Tested  by  Sutherland  and  Lingle 
(1972) 


Volume 

Young's  Modulus 

Poisson’s  Ratio 

Mass  Density 

Fraction  n(o) 

£fa) 

„<»> 

P™ 

^Tungsten 

0.022 

398  GPa 

0.28 

19194  kg/m3 

(2) 

'  'Aluminum 

0.978 

71.0  GPa 

0.34 

2700  kg/m3 
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Table  3.  Comparison  of  Effective  Moduli  of  a  Boron/ Aluminum  Unidirectionally  Fiber-Reinforced 
Composite  in  Units  of  10"  N/m 2 


Data3 

Mixture 

Model 

Square  Cell3 

Model 

Hexagonal  Cell3 

Model 

Eir' 

2.450 

2.551 

2.480 

2.551 

E}?> 

1.825 

1.868 

1.856 

1.872 

0.779 

0.661 

— 

0.661 

ElT} 

0.604 

0.578 

— 

0.578 

Ei «"* 

0.526 

0.604 

— 

0.606 

E^ 

0.566 

0.559 

0.451 

0.561 

*  After  Datta  and  f  edbeiter  /I9831 
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Figure  4.  Group  velocity  spectra  of  vavnref loctr.  ">cdo 
for  a  boron/epoxy  composite (Tauchert  and 
Guzelse, 1972) 
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Figure  5.  A  phase  velocity  spectrum  of  a  longitudinal 
wavereflect  mode  for  a  tungusten/aluminum 
composite  (Sutherland  and  Lingle,1972) 


